Small Ionized and Neutral Structures in the Diffuse Interstellar Medium 
ASP Conference Series, Vol. ** VOLUME**, 2007 
M. Haverkorn & W. M. Goss 



Pulsar Studies of Tiny-Scale Structure in the Neutral ISM 



J. M. Weisberg 

Department of Physics and Astronomy, Carleton College, Northfield, 
MN 55057 

S. Stanimirovic 

Radio Astronomy Lab, UC Berkeley, 601 Campbell Hall, Berkeley, CA 
94720 

Abstract. 

We describe the use of pulsars to study small-scale neutral structure in the 
interstellar medium (ISM). Because pulsars are high velocity objects, the pulsar- 
Earth line of sight sweeps rapidly across the ISM. Multiepoch measurements of 
pulsar interstellar spectral line spectra therefore probe ISM structures on AU 
scales. We review pulsar measurements of small scale structure in HI and OH 
and compare these results with those obtained through other techniques. 



1. Introduction 



The discovery of tiny-scale atomic structure (TSAS; lHeilesl ( 19971 )) in the neutral 



interstellar mediu m (ISM) was made in the VLBI neutral hydrogen (HI) absorp- 
tion experiment of Dieter et al. (1976)0] Meanwhile, a complementary technique 



for probing tiny-scale structure, using multiepoch pulsar observations, began in 
the mid-1980s. This technique takes advantage of pulsars' compact nature and 
their high velocities to probe the ISM on tiny scales by searching for temporal 
variations in spectral lines produced by the intervening medium. Until recently, 
these studies were confined to neutral hydrogen 21 cm absorption lines, but in 
the last few years the technique has been extended to investigations of OH in 
both absorption and stimulated emission. 



2. The Technique 

Pulsars subtend very small angles on the sky, even after their intrinsic size is 
broadened by interstellar scattering. They tend to be high-velocity objects, with 
transverse speeds of hundreds of km s _1 . Consequently the pulsar-Earth column 
probed by a pulsar signal is both needle thin and is dragged in short order across 
solar system size scales in the ISM. For example, a 100 km s _1 transverse pulsar 
velocity sweeps the line of sight across 20 AU in one year. Multiepoch pulsar 



1 See Brogan et al, this volume, for a review of interferometric studies of tiny-scale structure. 
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ISM measurements at reasonable intervals then probe ISM structure on these 
scales. 

In order to study the medium along the line of sight, one must measure the 
pulsar ISM spectrum^ at and near the frequency of a spectral line such as A ~ 21 
cm HI or A ~ 18 cm OH. The pulsar ISM spectrum, which is the spectrum of the 
pulsar alone (as modified by the intervening ISM) is not a direct observable since 
the telescope samples a complicated combination of signals throughout its beam. 
Some form of beam switching or spatial filtering is employed to separate the 
ISM spectrum of conventional steady sources. However, because of the pulsed 
nature of pulsar emission, it is straightforward in principle to derive the pulsar 
ISM spectrum via temporal switching; i.e., by differencing measurements of the 
spectrum when the pulsar is on, IpsR-On(v); and when it is off; IpsR—Offiy)'- 

IpSR-ISNl{v) = IpSR-Oniv) ~ IpSR-Off{ v )- (1) 

Ultimately we desire the fractional absorption or fractional stimulated emis- 
sion, IpsR-iSNiiy)/ 'Ipsr,o where Ipsr,o is the unabsorbed pulsar intensity de- 
termined from off-line spectral channels, or equivalently the (positive or nega- 
tive) optical depth t(v) of the gas along the pulsar-Earth line of sight, defined 
by 

t(v) = - hx(IpSR-ISM{v) / IpSR,o)- (2) 

Early efforts utilized one-bit spectrometers whose output was gated synchronously 
with the pulsar period into one of two buffers (pulsar-on or pulsar-q/f). Cali- 
bration of the spectra was difficult because o f limitations inherent in one-bit 
sampling of the wildly varying pulsar signal ( Weisberg et al.l [19801 ). Some of 



these limitations have been lessened through the years with the introduction of 
multibit spectrometers whose output can be directed into multiple pulsar phase 
bins rather than only two. The resulting multibit cube of data (intensity as 
a function of frequency and pulsar phase) can then be optimally processed to 
yield the pulsar-on, pulsar-q/f, and pulsar ISM spectra. Nevertheless, careful 
calibration remains crucial since the desired pulsar ISM spectrum represents the 
(usually small) difference between two very similar spectra. 



3. Neutral Hydrogen Measurements and Analyses 

As hydrogen is the most pervasive constituent of the ISM, it is natural to study 
fluctuations in HI. This section examines the pulsar and interferometric HI ex- 
periments and analyses. 

3.1. History of Pulsar Tiny-Scale HI Structure Experiments 

Early efforts with the pulsar HI ISM technique focused on measuring pulsar 
distances kinematically. This pr ocedure continues to ha v e great success, and 
provides primary calibration (e.g.. lFrail &: Weisberg! ( 1990l ); IWeisberg| ( 19961 )) for 



historically the relevant spectrum was called the "pulsar absorption" spectrum, but recent 
observations of pulsar-stimulated OH emission require that the term "absorption" be replaced 
when describing this technique generically. 
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mode ls of the galactic density distribution ( Taylor fc Cordeslll993l ; ICordes k, Lazio 
20021 ) . By the late 1980s, sufficiently accurate repeated measurements began to 



be made, and it was noticed that the pulsar ISM spectra changed over time in 
some cas es, suggesting c hange s in the tiny-scale structure of intevening gas. For 
example, IClifton et al.1 ([1988) found that the HI absorption spectrum of PSR 
B1821+05 changed significantly between 1981 and 1988, with the appearance 
at the l atter epoch of a p r eviou sly unobserved feature with r ~ 2 and Av ~ 1 



km s 1 . Deshpande et all ( 19921 ) showed that between ~ 1976 and 1981, HI ab- 



sorption toward B1154-52 did not change significantly; while toward B1557-50, 
a variation with At ~ 1 was interpreted as a cloud of size in the 1000 AU range. 

These early pulsar HI results inspired Frail, Weisberg, Cordes, &: Mathers 
(1994) to mount a dedicated multiepoch pulsar HI absorption experiment at 
Arecibo. Six pulsars were observed at three epochs, with time baselines ranging 
from 0.7 - 1.7 yr. These authors reported the presence of pervasive variations 
with At ~ 0.03 — 0.7, and associated HI column densities of 10 19 — 5 x 10 20 
cm -2 ; on scales of 5-100 AU. They indicated that 10-15% of cold HI is in the tiny 
structures, and additionally detected a correlation between absorption equivalent 
width variations and equivalent width itself. These results appeared to buttress 
the earlier VLBI findings of small scale structure, and they provided a strong 
impetus for further experimental and theoretical work. 

The theoretical work focused on explaining the very existence of TSAS. 
The presence of such structures is enigmatic, especially because t hey ap p ear to 
be significantly overpressured with respect to the ambient ISM. Heilesl ( 19971 ) 
emphasized the seriousness of the overpressure problem, and suggested that 
nonspherical cloud geometries such as filaments and sheets oriented w ith their 
long axes along the line of sight could ameliorate it. iDeshpandel ( 2000l ) argued 
that the observed changes do not result from discrete structures, but rather are 
merely the v ariations expec ted from a reasonable power-law spectrum of opacity 
fluctuations. Gwinn ( 200ll ) showed that velocity gradients in the HI, combined 
with interstellar scattering of the pulsar signals, could explain the variations. 



3.2. Recent Pulsar HI TSAS Measurements 



T he recent era of p ulsar TSAS experiments began with the Parkes observations 
of I Johnston et all ( 20031 ). Surprisingly, these investigators found no significant 
optical depth variations in their three-epoch, 2.5-yr observations of three pulsars, 
and were able to place an upper limit on column density variations toward 
PSR B1641-45 of 10 19 cm -2 , significantly below the Frail, Weisberg, Cordes, & 
Mathers (1994) detections. They showed that the earlier experiment did not 
fully account for the large increases in noise in absorption spectra at the line 
frequency, so that some of the apparently sign ificant variations actua lly were not. 
While no significant variations were seen by I Johnston et aL ( 20031 ) during the 
2.5 years of their TSAS experiment, they detected variations in the spectrum of 
PSR B1557-50 when compared with measurements ma de five years earlier. Thi s 
is the same pulsar whose spectrum was noted earlier bv lDeshpande et"ail ( 19921 ) 
to vary on similar timescales in the late 1970s. In comb i ning the results from 
four measurements over twenty-five years. Ijohnston et al.l ( 20031 ) concluded that 

the cloud causing the variations is ~ 1000 AU in size, with a density of ~ 10 4 
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Minter et al. ( 20051 ) performed a particularly exhaustive TSAS study on 
PSR B0329+54, which is very bright and almost circumpolar at the Green Bank 
Telescope. The investigators observed the pulsar continuously for up to 20 h in 
eighteen observing sessions over a period of 1.3 yr. They detected no HI opacity 
variations (At < 0.026 in most cases) for pulsar transverse offsets ranging from 
0.005 - 25 AU. 



The New Arecibo Pulsar HI TSAS Experiment Recently, we and our colleagues 
instigated a new set of Arecibo multiepoch pulsar H I TSAS studies. A dvances 
in digital technology since the earli er observations oflFrail et al" ( 19941 ) allowed 
the use of a multibit spectrometer ( Jenet et al.ll200ll ) whose calibration should 
minimize systematic errors; and re ceiver improvements decrea sed the noise. Pre- 
liminary results were published by Stanimirovic et al.1 (|2003al ) and a final report 
is near ing complet i on dStanimirovic et al.ll2007l ). We observed the same six pul- 
sars as iFrail et all (|l99i T at four epochs, yielding time baselines between 0.2 and 
1.3 yr, with pulsars moving between 1 and 200 AU between sessions. The line 
of sight parameters for the six pulsars are given in Table 1. 




Figure 1. Differences between p airs of HI absorption spect ra for B1929+10. 
Data were gathered at Arecibo by IStanimirovic et al.l (|2007f ). The horizontal 
bar delineates the velocity range where significant absorption is present in 
single-epoch spectra. Dashed lines show ±2<r error envelopes. Note that only 
a few features extend outside the envelopes. 
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Table 1. Line of sight parameters in the Arecibo pulsar HI TSAS study 
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(deg) 


b 
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^transverse 

(AU / yr) 
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3.5 
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B0823+26 
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0.4 


40 
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69 
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130 
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B1737+13 


37 


22 


4.7 
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J1932+1059 


B1929+10 


47 


-3.9 


0.3 


40 


J2018+2839 


B2016+28 


G8 


-4.0 


1.0 
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For each pulsar, we differenced the six pairs of absorption spectra corre- 
sponding to the four observing sessions. We found few statistically significant 
optical depth variations. The rare cases of significant changes were mostly seen 
toward PSR B1929+10 and PSR B2016+28 (see Figs. 1 and 2). The detections 
result from variations having column densities of ~ 10 8—19 cm -2 , with implied 
volume densities (for spherical clouds) of ~ 10 4-5 cm" 3 . It is worth noting that 
the PSR BT929+10 fluctuations have a factor of 3-4 lower column densities than 
the "canonical" HI TSAS discussed bv lHeiled (|l997l b 

We have also studied the integrated properties of the variations using equiv- 
alent width EW, where EW= / rdv. Fig. 3 shows |AEW| versus EW for all 
our pulsars except B0540+23, whose optical depth integrals become very uncer- 
tain near r ~ oo features. PSR 1929+10 exhibits the only significant (> 3a) 
nonzero values of |AEW|. Our results on PSRs B1737+13 and particularly 
B2016+28 [primarily upper limits with an occasional marginal (< 2a) detection] 
are far lower than the detections of Frail, Weisberg, Cordes, & Mathers (1994) 
[|AEW| < 0.2 vs. |AEW| ~ 0.5 for B1737+13 and |AEW| < 0.3 vs. |AEW| ~ 
(1 — 5) for B2016+28]. We believe that at least some of the variations they saw 
were due to small calibration errors. Our results do not rule out their finding of 
a correlation between |AEW| and EW, but the evidence for such a trend appears 
marginal in our data. 



3.3. A Sketch of Recent Interferometric HI TSAS Experiments 



The state of the art of interferometric measure ments of HI TSAS has also ad- 
vanced significantly since the first detection of iDieter et al.l ( 19761 ). especially 



with the advent of the VLBA. The sou r ces 3 C 147 and 3C 138 show the most 
prominent fluctuations. iFaison Gosd ( 200ll ). who have done the most recent 



observations and analyses of 3C 147, find Ar_~ 0.3 over ~ 25 AU scales. The 



newest results on 3C 138 bv lBrogan et al.l ( 20051 ) indicate variations of At ~ 0.5 
on 25 AU scales. The fluctuations appear in approximately 10% of their pixels 
on the source. Most modern interferomete r results on other s ources yield smaller 
and less frequent optical depth variations ( Faison et al.|[i~998l ). See the contribu- 
tion by Brogan et al in this volume for a thorough review of these measurements. 
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Figure 2. Differences between pairs of absorption spectra for B2016+28. 
See Fig. 1 for details. 



3.4. Discussion of HI TSAS Results 

The various pulsar and interferometric measurements reviewed above appear 
to lead to divergent conclusions about the properties and pervasiveness of HI 
TSAS. Yet Occam's Razor suggests that we search for the simplest explanation 

that can explain them all. 

Among the pulsar measurements, it appears that iFrail et aL I (119941 ) . who 



found that TSAS is rather pervasive, provide the principal outliers. We have 
noted above that there appear to be calibration issues in these measurements 
which may have led to overestimates of the prevalence of optical depth varia- 
tions. In our subsequent analyses, we have replaced their res ults with our newer 
mea surements of the sam e pulsars ( Stanimirovic et al. 20071 ). 



Brog an et al. ( 20051 ) note that the two sources exhibiting the most signif- 



icant evidence for TSAS among HI interferometer results, 3C 138 and 3C 147, 
are also the only two with angular sizes significantly larger than the detected 
gas fluctuations. They propose that the lower measured levels of variations in 
the other cases results from a selection effect due to incomplete sampling of the 
relevant angular and spatial scales. They also attempt to model pulsar results by 
studying the statistics of one-dimensional spatial cuts across the face of 3C138. 
They conclude that the relatively low level of recent pulsar TSAS detections 
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Figure 3. Equivalent width changes |AEW| versus EW in the pulsar HI 
experiment of Stanimi rovic et all (|2007T) . where equivalent width EW= / rdv. 



can also be explained by an incomplete sampling of angular scales due to the 
episodic nature of the observations. However, their model pulsar observations 
are all separated by the same angular scale, namely the characteristic scale of 
the HI fluctuations; whereas multiepoch observations of a given pulsa r actu - 
ally sample a varie ty of scales. The rec ent observations of iMinter et al.l ((2005), 
published after the lBrogan et al.l ( 2005I ) analysis was complete, are particularly 
notable in finding no fluctuations over a very wide and almost continuous range 
of scales. 

In an attempt to make progress toward a comprehensive explanation of 
the various observations, we have plotted HI TSAS detections and upper limits 
as a function of their spatial scale in Fig. 4. The expect ed level of optica l 
depth variations as a function of scale calculated from the iDeshpandel ( 20001 ) 
power law fluctuation model is also shown as a sloping line. As these values 
are extrapolated from much larger observed scales, they could be significantly 
in error but the slope is more secure. 

We note that the iBrogan et al. ( 20051 ) interferometric detection is at a 
higher level than any other, which suggests that the line of sight may indeed 
have a higher level of fluctuations than is typical. Note however that several of 
our pulsar detections at slightly smaller scales (mostly PSR 2016 +28), though 



exhibi ting somewhat l ower optical depth s, are consistent with the lBrogan et al 



(|2005T ) detection if the iDeshpandl (|2OO0h power law slope is correct. Meanwhile 
numerous other pulsar detections and upper limits lie well below th e above men- 
tioned ones; most especially our PSR B 1929+ 10 results and the IMinter et al 



(2005) upper limits over a wide range of scales. The figure also emphasizes that 
pulsar TSAS experiments have made detections on many spatial scales, whereas 
the two most secure interferometric results both yield a single scale of ~ 25 AU 
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Figu re 4. HI op t ical d epth fluctuations as a function of spatial scale. 
The iBrogan et all (l2005f ) interferometric 3C 138 fluctuations (star), the 
Uohnston et al.1 (120031) pulsar det ection (filled square) and upper limits (open 
squares), the Min ter et al.1 (120051) p ulsar upper limits (dot-dashed horizontal 
line), and IStanimirovic et al.1 (|2007l ) pulsar det ections (filled circ les) and up- 
per limits (crosses) are shown, along with the iDeshpandd (|2000| ) theoretical 
expectation (sloping line) . The depth of the fluctuation in units of <j is placed 
next to each detection. 



( Broean et al.1 120051 ). Hence we conclude that the observed variations in level 



and scale of HI TSAS are more than a selection effect, and that they are strongly 
dependent up on local conditions along a par ticular line of sight. 

Recently iBraun &: Kanekarl (2005|) and IStanimirovic fc Heilesl (|2005l ) have 
detected a population of small (10 3-4 AU) HI clouds via very sensitive absorp- 
tion measurements of extragalactic background sources. The optical depths and 
column densities range below 1 % and 10 18 cm -2 , respectively. It is not yet clear 
whether these features represent a distinct population or merely an extension of 
the pulsar-detected TSAS toward larger spatial size and smaller optical depth. 

We have examined the lines of sight toward our two pulsars with the most 
significant fluctuations. Nothing particular was found in the direction of PSR 
B2016-I-28. However, the line of sight toward PSR B1929+10 appears to run 
al ong the wall of the Loca l Bubble for a significant stretch, according to the map 
of Lallement et al.l ( 20031 ) . While the Local Bubble itself is probably a relatively 
quiescent region, its bo undary is thought to be marked by enhanced turbulence 
( Phillips Sz Clega 19921 ). We currently regard this finding as potentially impor- 
tant but tentative, and we await further refinements in knowledge of the location 
of the Local Bubble boundary. 
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4. Studies of Molecular Tiny-Scale Structure 

The tiny-scale structure of molecular gas in the ISM is an interesting area for 
investigation, both in its own right and in comparison with neutral atomic hydro- 
gen (§3). Interstellar molecular spectral lines usually have much smaller optical 
depths than HI, making them harder to detect. Therefore most successful radio 
observations of small scale structure in molecular gas awaited relatively recent 
adv ances in sensitiv i ty and ang ular resolution. 

iMarscher et~aH (|l993l ) and lMoore & Marscheil (|l995l ) did a series of molec- 



ular absorption line experiments using extragalactic background sources. They 
detected significant time-variable 4.8 GHz H2CO absorption with At = 0.02 — 
0.03 in front of two compact sources, which they interpret as resulting from the 
superposition of several clumps along the line of sight whose size is ~ 10 AU or 
less, having HI column densities of ~ 10 20 cm -2 and volume densities of ~ 10 6 
cm~ 3 . 

Pulsars can be excellent probes of molecular gas for all the reasons given 
in §2, providing only that the pulsar signal is sufficiently strong at the line 
frequency. However, pulsars have rather steep spectra, so the only lines matching 

this criterion are the four 18 cm lines of OH. 

The first p ublished pulsar OH ISM experiments were done bv lSlvshl (|l972l ) 



and lGaltl ( 1974! ). Both of these investigators detected no absorption in the spec- 



trum of PSR B0329+54; while Slysh detected no absorption in several others. 
His detection of OH absorption at 1667 MHz in PSR B1749-28 was not confirmed 
by later observations. 

Taking advantage of the high sensitivity of Are cibo, we and our colleague s 
did the first modern pulsar studies of interstellar OH ( Stanimirovic et al. 2003b). 



Among the seven pulsars studied, only PSR B1849+00 clearly exhibits OH ab- 
sorption. This line of sight is remarkable in that it passes near the edge of the 
SNR Kes 79 (G33. 6+0.1). Fig. 5 displays our OH spectra gathered in this 
direction. The absorption lines in the pulsar-off spectra (bottom row) result 
from SNR continuum radiation being absorbed by foreground gas. Note that 
these pulsar-off absorption lines exhibit much smaller optical depths than do the 
pulsar OH spectra (top row) at v ~ 100 km s _1 (r ~ 0.05 vs. 0.4 at 1665 MHz 
and r ~ 0.02 vs. 0.9 at 1667 MHz). These smaller pulsar-off optical depths are 
much more representative of typical OH absorption optical depths in the plane 
(e.g., Dickey et al 1981), so that our measur ed pulsar OH optical depths are 
anomalously high. IStanimirovic et ail ( 2003bl ) discuss several possible explana- 



tions for the high optical depths in the pulsar OH interstellar spectra. The most 
likely is that the tiny column through the ISM sampled by the pulsar signal 
has encountered a small (< 15"), high density (> 10 5 cm -3 ) molecular cloudlet. 
Meanwhile, the pulsar-off spectrum, in sampling absorption across the full tele- 
scope beam, represents a solid-angle average across a clumpy medium consisting 
of relatively high optical depth, small molecular cloudlets, embedded in a lower 

den sity medium. 

IWeisberg et al.l ( 20051 ) continued the search for additional pulsar OH in- 



terstellar lines with the Parkes telescope. Of eighteen pulsars searched, only 
B1641-45 exhibited detectable pulsar ISM spectra (see Fig. 6). Note that the 
pulsar ISM spectra (right column) exhibit significantly larger optical depths 
than do the corresponding pulsar-off spectra (left column), continuing the trend 
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Figure 5. Arec i bo O H spectra toward PSR E-1849+00, after 
Stan imirovic et aT] ((2003b). The top row shows the pulsar ISM spec- 
tra at 1665 and 1667 MHz. The bottom row, taken when the pulsar is off, 
shows emission and absorption throughout the telescope beam, which is filled 
by the SNR Kes 79 (G33.6+0.1). Compare pulsar ISM optical depths (top) 
with pulsar-off, SNR optical depths (bottom). 



discussed above in PSR B1849+00. So it again appears that the pulsar-Earth 
pencil beam has intercepted dense cloudlets in the pulsar ISM spectra, whose 
properties suffer angular dilution in the pulsar-off spectra. Note also that stimu- 
lated emission is present in the 1720 MHz pulsar ISM spectrum. This is the first 
detection of a pulsed interstellar maser. Since the stimulated emission appears 
in the spectrum of the pulsar alone, it represents a direct demonstration that 
the emission is sti mulated by pulsar photons. 

iMinterl (|2005l ) has also recently searched for OH lines in the spectra of five 
pulsars with the GBT. He detected absorption toward PSR B1718-35. 



4.1. Discussion of Molecular Results 



Theory suggests that mole c ular clouds are a dy n amic, turbulent e nviron ment 



(|Mac Low fc Klessenl 120041 : iBonnell et all l200fil i iBrunt fc HeveH (120021 ) per 



formed a principal component analysis of CO maps, finding a slope of 2.17 for 
the energy spectrum, which they interpreted as evidence for continuous energy 
injection at small scales. Thus it is not surprising that both the formaldehyde 
results and our groups' pulsar OH work indicate the presence of tiny-scale struc- 
ture in molecular gas. Our conclusions result from our detection of larger optical 
depths in the pulsar ISM spectra than in the pulsar-off spectra. Multiepoch ex- 
periments are now in progress to search for time variability in the these lines, 
much like the pulsar HI experiments of §3.1 and 3.2. Detection of time variability 
will provide powerful additional support for the existence of TSAS in OH. 
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PSR-OFF, 1720 MHz 




PULSAR, 




1720 
MHz 






Figure 6. Pulsar-off (left) and Pulsar ISM (right) spectra a t the four 18 cm 
OH lin es toward PSR B1641-45, from Parkes observations bv lWeisberg et ahl 
(2005). All spectra are plotted with the same optical depth scale. Compare 
pulsar-off and pulsar ISM optical depths, and also note the pulsed interstellar 
OH maser stimulated by the pulsar at 1720 MHz. 
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